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Abstract 

We present a preliminary direct measurement of the parity violating coupling of the Z° 
to strange quarks, A B , derived from a sample of approximately 300,000 hadronic decays of 
Z° bosons produced with a polarized electron beam and recorded by the SLD experiment 
at SLAC between 1993 and 1997. Z° — > ss events are tagged by the presence in each 
event hemisphere of a high-momentum if , K® or A°/A° identified using the Cherenkov 
Ring Imaging Detector and/or a mass tag. The CCD vertex detector is used to suppress 
the background from heavy flavor events. The strangeness of the tagged particle is used 
to sign the event thrust axis in the direction of the initial s quark. The coupling A s 
is obtained directly from a measurement of the left-right-forward-backward production 
asymmetry in polar angle of the tagged s quark. To reduce the model dependence of the 
measurement, the background from uu and dd events is measured from the data, as is the 
analyzing power of the method for ss events. We measure: 

A s = 0.82 ± 0.W(stat.) ± 0. 08 (syst.) (preliminary). 
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1. Introduction 



Measurements of the fermion production asymmetries in the process e + e~ — > Z° — > ff 
provide information on the extent of parity violation in the coupling of the Z° boson 
to fermions of type /. At Born level, the differential production cross section can be 
expressed in terms of x = cos 6, where 9 is the polar angle of the final state fermion / 
with respect to the electron beam direction: 



a 



rlfjf 

ff x) = — oc (1 - A e P e )(l + x 2 ) + 2A f (A e - P e )x, 



where P e is the longitudinal polarization of the electron beam, the positron beam is 
assumed unpolarized, and the asymmetry parameters Af = 2vfdf/(Vf + aj) are defined 
in terms of the vector (vf) and axial- vector (a/) couplings of the Z° to fermion /. If 
one measures the polar angle distribution for a given final state //, one can derive the 
forward-backward production asymmetry, Ap B , which depends on both the initial and 
final state asymmetry parameters as well as on the beam polarization: 

_ A e -P e x 
Afb{X) ~ af(x) + af(-x) ~ ZAf l- A e P e 1 + x 2 ' { ) 

For zero polarization, one measures the product of couplings A e Aj. If one measures 
the distributions in samples taken with negative (L) and positive (R) beam polarization 
of magnitude P e , then one can derive the left-right-forward-backward asymmetry, A FB , 
which is insensitive to the initial state coupling: 

AU*) = ^ + ^~ X)) ~ + = 2\PMf-^. (3) 

(a f L (x) + ai(-x)) + + a f L (-x)) f I + 

A number of previous measurements have been made by experiments at LEP and 
SLC of A e , A^, A T , A c and At, In contrast, very few measurements exist for the light 
flavor quarks, due to the difficulty of tagging specific light flavors. It has recently been 
demonstrated experimentally that light flavored jets can be tagged by the presence 
of a high- momentum 'leading' identified particle that has a valence quark of the desired 
flavor, for example a K~ (K + ) meson could tag an s (s) jet. However the background 
from other light flavors (a u jet can also produce a leading K~), decays of B and D 
hadrons, and nonleading kaons in events of all flavors is large, and neither the signal nor 
the background has been well measured experimentally. 

The DELPHI collaboration has measured @ A S FB = 0.131 ±0.03h(stat.) ±0.013(syst.) 
and Ap S B = 0. 1 12 ±0. 031 (stat. ) ±0.054(syst.), respectively. However the extraction of the 
asymmetry parameters from the measured production asymmetries is model dependent. 
The OPAL collaboration has measured || A U FB = 0.040 ± 0.0Q7(stat.) ± 0.028(syst.) and 
Ap S B = 0.068 ± 0.035(stat.) ± 0. 011 (syst.). In their measurement, most of the background 
contributions and analyzing powers are determined from double-tagged events in the data. 
This eliminates most of the model dependence, but results in limited statistical precision. 
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We present a direct measurement of the asymmetry parameter for strange quarks, 
A s , using a sample of 300,000 hadronic Z° decays recorded by the SLD experiment at the 
SLAC Linear Collider between 1993 and 1997, with an average electron beam polarization 
of 74%. We use identified strange particles to tag s and s jets. The analyzing power 
of the tags for true ss events, as well as the relative contribution of uu + dd events 
were determined from the data. This procedure removes much of the model dependence, 
yielding an error that is statistically dominated. 



2. Hadronic Event and ss Event Selection 

A general description of the SLD can be found elsewhere 0. The trigger and initial 
selection criteria for hadronic Z° decays are described in Ref. ||. In order to reduce 
the effects of decays of heavy hadrons, we selected light flavor events (uu, dd and ss) by 
requiring each high-quality [0] track in the event to have a transverse impact parameter 
with respect to the IP of less than three times its estimated error. Finally, the CRID 
(Cherenkov Ring Imaging Detector) was required to be operational. The selected sample 
comprised roughly 94,000 events with an estimated background contribution of 11% from 
cc events, 2% from bb events, and a non-hadronic background contribution of 0.10±0.05%, 
dominated by Z° — > t + t~ events. 

For the purpose of estimating the efficiency and purity of the event flavor tagging and 
the particle identification, we made use of a detailed Monte Carlo (MC) simulation of the 
detector. The JETSET 7.4 || event generator was used, with parameter values tuned 
to hadronic e + e~ annihilation data |§, combined with a simulation of S-hadron decays 
tuned |TIJ to T(4S) data and a simulation of the SLD based on GEANT 3.21 [ HJ . 



Then high-momentum strange particles are selected. The CRID allows K to be 
separated from p/p and ti^ with high purity for tracks with p > 9 GeV/c as described 
in detail in [12|]. A track is tagged as a K if the log-likelihood for this hypothesis 
exceeds both, the ^ and the p/p log-likelihoods, by at least 3 units. The average purity 
of the sample was estimated using the simulation to be 89%. K® and A°/A° are 
reconstructed |T2| in the modes K® — > 7r + 7r~ and A°(A°) — > p(p)n T and are identified 



by their long flight distance, reconstructed mass, and accuracy of pointing back to the 
primary interaction point, and are required to have p > 5 GeV/c. Pairs of tracks with 
invariant mass m nn within 12 MeV/c 2 of the nominal K® mass are identified as K®. The 
simulation predicts that the average purity of the sample is 91%. In the case of the 
A°/A°, information from the Cherenkov Ring Imaging Detector is used to identify the p/p 
candidate. A°/A° are identified by requiring the invariant mass of pairs of tracks, m pn , 
to be within 5 MeV/c 2 of the nominal A°/A° mass and a combination of K® rejection 
and particle identification for the high-momentum track. The simulation predicts that 
the average purity of the A°/A° sample is 84%. 

These strange particles are then used to tag s and s jets as follows. The event is 
divided into two hemispheres by a plane perpendicular to the thrust axis. We require 
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each of the two hemispheres to contain at least one identified strange particle (K , K® 
or A°/A°); for hemispheres with multiple strange particles we only consider the one with 
the highest momentum. We require at least one of the two hemispheres to have definite 
strangeness (i.e. to contain a K ± or A°/A°). In events with two hemispheres of definite 
strangeness, the two hemispheres are required to have opposite strangeness (e.g. K + K~). 
This procedure increases the ss purity substantially compared with a single tag; thus, for 
these events, the model dependence of the measurement is reduced. 

Table [I] summarizes the composition of the selected event sample for data and simula- 
tion for each of the 5 tagging modes used. The number of events for each mode shown is 
in good agreement with the MC prediction. The ss purity and ss analyzing power were 
estimated from the simulation. 



Table 1: Summary of selected event sample for 5 modes in data and simulation. 



Mode 


Data 


MC prediction 


ss purity 


ss analyzing power 


K+K- 


619 


620 


0.76 


0.94 




86 


82 


0.65 


0.86 


A U A U 


1 


6 


0.59 






502 


531 


0.64 


0.68 


A»K U S ,A°K" S 


52 


52 


0.54 


0.44 


Total: 


1260 


1291 


0.69 


0.82 



The combined ss purity of all modes is 69%, and the predicted background in the 
selected event sample consists of 10% uu, 9% dd, 11% cc, and 1% bb events. The analyzing 
power is defined as a s = ° right ~ L rong where N r s l9tl (N™ ron9 ) denotes the number of ss 

jV s -j-iV s 

events in which a particle of negative strangeness is found in the true s(s) hemisphere. 
The average analyzing power for all modes is predicted by the simulation to be 0.82. The 
K ± K T mode has a substantially higher analyzing power and ss purity than the other 
modes. The initial s quark direction is approximated by the thrust axis, t of the event, 
signed to point in the direction of negative strangeness, x = cos9 s = S-^i z , where S and 
p denote the strangeness and the momentum of the tagging particle. 

Figure [I] shows the polar angle distributions, for all 5 tagging modes combined, of the 
signed thrust axis for left-handed (P e < 0) and right-handed (P e > 0) electron beams. 
The expected production asymmetries, of opposite sign for the left-handed and the right- 
handed beams, are clearly visible. 
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Figure 1: Polar angle distributions of the thrust axis, signed to point in the direction of 
negative strangeness, of the tagged strange particle, for negative (left) and positive (right) 
beam polarization. The dots show data and the histogram shows our fit to the data. Our 
estimates of the non-ss backgrounds are indicated by the hatched histograms. 

3. Extraction of A s and Systematic Uncertainties 

A s is extracted from these distributions by a binned maximum likelihood fit. The fitting 
function is given by: 

P(x) = D(x) N f (l + x 2 + 2(1 + 5)a f A f A z x). (4) 

The function D(x) describes the acceptance and the strange particle identification ef- 
ficiencies. Nf = N events Rfef denotes the number of events in the sample of flavor / 
(/ = u, d, s, c, b) in terms of the number of selected hadronic events N events , Rf = T(Z° — > 
ff)/T(Z° — > hadrons) and the tagging efficiencies e/; 5 = —0.013 corrects for the effects of 
hard gluon radiation ||13|| ; a/ denotes the analyzing power for tagging the / rather than the 
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/ direction; Af is the asymmetry parameter for flavor /; and Az = {A e — P e ) /(l — A e P e ). 
The function D(x) was calculated from the simulation and verified by comparing data 
and simulated x-distributions of identified K ± , Kg, and A /A . The parameters e c , e&, 



and a c , ay, for the heavy flavors are taken from the MC simulation . The world average 



experimental measurements of the parameters A c , At, R c , Rb [|1| were used. For the light 
flavors, the relevant parameters in the fitting function are derived where possible from 
the data as described below. The total number of light flavor events, N U( i s , is determined 
by subtracting the simulated number of heavy flavor events from the entire event sample. 
The ratio N u d/N s is 0.27 ± 0.03. The asymmetry parameters A u and Ad are set to the 
Standard Model values. The ss analyzing power, a s , is 0.82 ± 0.03 (averaged over all 
modes). The combined (uu + dd) analyzing power, a u a, is estimated to be —0.41 ± 0.24. 
The result of the fit is shown as a histogram in Figure [1]. The fit quality is good with 
a x 2 °f 12.9 for 24 bins. Also included are our estimates of non-ss background. The 
cross-hatched histograms indicate cc + bb backgrounds which are seen to show asymme- 
tries of the same sign and similar slope to the total distribution. The hatched histograms 
indicate uu + dd backgrounds showing asymmetries of the opposite sign and slope to the 
total distribution. The A s value extracted from the fit is A s = 0.82 ± 0A0(stat.). 

The understanding of the parameters used as inputs to the fitting function and of their 
uncertainties is crucial to this analysis. The characteristics of heavy flavor events relevant 
to this analysis have been measured experimentally, and our simulation |j, ||], [Tl| has 
been tuned jl0] to reproduce these results. Standard systematic variations of the heavy 
flavor simulation were considered, and the resulting uncertainties are a small contribution 
to the total systematic error. Other small contributions to the systematic error include 
those from the 0.6% uncertainty in the correction for the effect of hard gluon radiation, 
and the 0.8% uncertainty in the beam polarization. For the light flavors, there are few 
experimental constraints on the relevant input parameters. Qualitative features such as 
leading particle production || , short range rapidity correlations between high-momentum 
KK and baryon-antibaryon pairs [[14]] and long-range correlations between several particle 
species 0] have been observed experimentally, but these results are not sufficient to 
quantify the analyzing power of the strange-particle tag or the uu and dd background. 
Our MC simulation provides a reasonable description of the above observations and was 
used to evaluate the parameters used in the fit, but was calibrated using data. For the 
analyzing power in ss events, we investigated the rate of production of wrong-sign kaons 
by counting events in which we find three identified K ± and/or K®. Since we found that 
the MC prediction for the number of 3-kaon jets is consistent with the data, we used the 
simulated a s = 0.82 as our central value for the analyzing power in ss events and the 
statistical error on the data-MC comparison as uncertainty. We also counted hemispheres 
containing a K + K + or K~K~ pair, obtaining consistent but less precise constraints. 

For the calibration of the relative uu + dd background level, N u d/N s , we counted the 
number of hemispheres containing an identified K + -K~~ pair or an identified K^-K pair. 
The MC prediction is consistent with the data and was used as the central value. Another, 
less precise contraint was obtained from events that were tagged by kaons of the same 
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Table 2: Summary of systematic uncertainties. 



Source 


Comments 


Systematic variation 


8A S /A S 


heavy flavor modelling 


MC/world averages 


Ref. HI, EL [TOll 


0.008 


hard gluon radiation 


Stav-Olsen with 
bias correction 


(1.3 ±0.6)% 


0.006 


beam polarization 


data 


(73.7 ±0.8)% 


0.008 


a s 


MC constrained by 
3-K jets in data 


0.82 ±0.03 


0.038 


O-ud 


-a s < a ud < 


-0.41 ±0.24 


0.071 


A U( i 


Standard Model 






N ud /N s 


MC constrained by 
2-K jets in data 


0.27 ±0.03 


0.046 


MC statistics 






0.021 


Total: 






0.096 



sign in both hemispheres. The above checks are also sensitive to the analyzing power of 
uu ± dd events, a ud . However, with the present event statistics we cannot obtain a tight 
constraint on this quantity. We therefore assume that a u d must be negative, since u and 
d jets must produce a leading K + rather than K~ , and that \a U( i\ must be less than a s , 
since there is always a companion particle of opposite strangeness in a u or d jet that 
will tend to dilute the analyzing power. We take these as hard limits, —0.82 < a u d < 0, 
use the middle of the range for our central value and assign an uncertainty equal to the 
range divided by \/l2. The simulation predicts a value of a u d = —0.38, consistent with 
our estimate. Table ||| summarizes the systematic uncertainties. The total systematic 
uncertainty is SA S = 0.08. 

4. Summary and Conclusion 

We have presented a preliminary direct measurement of the parity violating coupling of 
the Z° to strange quarks, A s , derived from a sample of approximately 300,000 hadronic 
decays of Z° bosons produced with a polarized electron beam and recorded by the SLD 
experiment at SLAC between 1993 and 1997. The coupling A s is obtained directly from a 
measurement of the left-right-forward-backward production asymmetry in polar angle of 
the tagged s quark. The background from uu and dd events is measured from the data, 
as is the analyzing power of the method for ss events. A binned maximum likelihood fit 
is used to obtain the result: 

A s = 0.82 ± 0A0(stat.) ± 0.08(syst.) (preliminary). (5) 
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This result is consistent with the Standard Model expectation for A s . Our measurement 
can be used to test the universality of the coupling constants by comparing it with the 
world average value for ffl. The two measurements are consistent. 

In order to compare with previous measurements of A S FB and A d p S B (see section 1), we 
must assume a value of A e . Using A e = 0.1499 [|1|] and neglecting the small uncertainty 
on A e , the DELPHI measurements translate into A s = 1.165 ± 0.311(stat.) ± 0. 116 (syst.) 
and A d jS = 0.996 ± 0.276(stat.) ± 0A80(syst.). Similarly, the OPAL measurement yields 
Ad,s — 0.605 ± 0.311(stat.) ± 0.098(syst.). Our measurement is consistent with these and 
represents a substantial improvement in precision. 
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